Drug-resistance profiling and transmission dynamics of multidrug-resistant Mycobacterium tuberculosis in Saudi Arabia revealed by whole genome sequencing.

BACKGROUND: In Saudi Arabia, cross-border transmission of multidrug-resistant (MDR) Mycobacterium tuberculosis complex (MTBC) strains might be particularly fostered by high immigration rates. Herein, we aimed to elucidate the transmission dynamics of MDR-MTBC strains and reveal a detailed prediction of all resistance-conferring mutations for the first- and second-line drugs. METHODS: We investigated all MDR-MTBC strains collected between 2015 and 2017 from provincial mycobacteria referral laboratories and compared demographic and clinical parameters to a cohort of non-MDR-TB patients using a whole genome sequencing approach. Clusters were defined based on a maximum strain-to-strain genetic distance of five single-nucleotide polymorphisms (SNPs) as surrogate marker for recent transmission, and then investigated molecular drug-resistance markers (37 genes). RESULTS: Forty-eight (67.6%) MDR-MTBC strains were grouped in 14 different clusters, ranging in size from two to six strains; 22.5% (16/71) of all MDR-MTBC isolates were predicted to be fully resistant to all five first-line drugs, and five strains (7.0%) exhibited fluoroquinolone resistance. Moreover, we revealed the presence of 12 compensatory mutations as well as 26 non-synonymous SNPs in the rpoC gene and non-hotspot region in rpoB, respectively. CONCLUSION: Optimized TB molecular surveillance, diagnosis, and patient management are urgently needed to contain MDR-MTBC transmission and prevent the development of additional drug resistance.

Characterizing the morbid genome of ciliopathies.

BACKGROUND: Ciliopathies are clinically diverse disorders of the primary cilium. Remarkable progress has been made in understanding the molecular basis of these genetically heterogeneous conditions; however, our knowledge of their morbid genome, pleiotropy, and variable expressivity remains incomplete. RESULTS: We applied genomic approaches on a large patient cohort of 371 affected individuals from 265 families, with phenotypes that span the entire ciliopathy spectrum. Likely causal mutations in previously described ciliopathy genes were identified in 85% (225/265) of the families, adding 32 novel alleles. Consistent with a fully penetrant model for these genes, we found no significant difference in their "mutation load" beyond the causal variants between our ciliopathy cohort and a control non-ciliopathy cohort. Genomic analysis of our cohort further identified mutations in a novel morbid gene TXNDC15, encoding a thiol isomerase, based on independent loss of function mutations in individuals with a consistent ciliopathy phenotype (Meckel-Gruber syndrome) and a functional effect of its deficiency on ciliary signaling. Our study also highlighted seven novel candidate genes (TRAPPC3, EXOC3L2, FAM98C, C17orf61, LRRCC1, NEK4, and CELSR2) some of which have established links to ciliogenesis. Finally, we show that the morbid genome of ciliopathies encompasses many founder mutations, the combined carrier frequency of which accounts for a high disease burden in the study population. CONCLUSIONS: Our study increases our understanding of the morbid genome of ciliopathies. We also provide the strongest evidence, to date, in support of the classical Mendelian inheritance of Bardet-Biedl syndrome and other ciliopathies.

Characterization of Greater Middle Eastern genetic variation for enhanced disease gene discovery.

The Greater Middle East (GME) has been a central hub of human migration and population admixture. The tradition of consanguinity, variably practiced in the Persian Gulf region, North Africa, and Central Asia, has resulted in an elevated burden of recessive disease. Here we generated a whole-exome GME variome from 1,111 unrelated subjects. We detected substantial diversity and admixture in continental and subregional populations, corresponding to several ancient founder populations with little evidence of bottlenecks. Measured consanguinity rates were an order of magnitude above those in other sampled populations, and the GME population exhibited an increased burden of runs of homozygosity (ROHs) but showed no evidence for reduced burden of deleterious variation due to classically theorized 'genetic purging'. Applying this database to unsolved recessive conditions in the GME population reduced the number of potential disease-causing variants by four- to sevenfold. These results show variegated genetic architecture in GME populations and support future human genetic discoveries in Mendelian and population genetics.

Revisiting disease genes based on whole-exome sequencing in consanguineous populations.

Assigning a causal role for genes in disease states is one of the most significant medical applications of human genetics research. The requirement for at least two different pathogenic alleles in the same gene in individuals with a similar phenotype to assign a causal link has not always been fully adhered to, and we now know that even two alleles may not necessarily constitute sufficient evidence. Autozygosity is a rich source of natural "knockout" events by virtue of rendering ancestral loss-of-function (LOF) variants homozygous. In this study, we exploit this phenomenon by examining 523 exomes enriched for autozygosity to call into question previously published disease links for several genes based on the identification of confirmed homozygous LOF variants in the absence of the purported diseases. This study highlights an additional advantage of consanguineous populations in the quest to improve the medical annotation of the human genome.

Autozygome sequencing expands the horizon of human knockout research and provides novel insights into human phenotypic variation.

The use of autozygosity as a mapping tool in the search for autosomal recessive disease genes is well established. We hypothesized that autozygosity not only unmasks the recessiveness of disease causing variants, but can also reveal natural knockouts of genes with less obvious phenotypic consequences. To test this hypothesis, we exome sequenced 77 well phenotyped individuals born to first cousin parents in search of genes that are biallelically inactivated. Using a very conservative estimate, we show that each of these individuals carries biallelic inactivation of 22.8 genes on average. For many of the 169 genes that appear to be biallelically inactivated, available data support involvement in modulating metabolism, immunity, perception, external appearance and other phenotypic aspects, and appear therefore to contribute to human phenotypic variation. Other genes with biallelic inactivation may contribute in yet unknown mechanisms or may be on their way to conversion into pseudogenes due to true recent dispensability. We conclude that sequencing the autozygome is an efficient way to map the contribution of genes to human phenotypic variation that goes beyond the classical definition of disease.

Identification of the tetraspanin CD82 as a new barrier to xenotransplantation.

Significant immunological obstacles are to be negotiated before xenotransplantation becomes a clinical reality. An initial rejection of transplanted vascularized xenograft is attributed to Galα1,3Galß1,4GlcNAc-R (Galα1,3-Gal)-dependent and -independent mechanisms. Hitherto, no receptor molecule has been identified that could account for Galα1,3-Gal-independent rejection. In this study, we identify the tetraspanin CD82 as a receptor molecule for the Galα1,3-Gal-independent mechanism. We demonstrate that, in contrast to human undifferentiated myeloid cell lines, differentiated cell lines are capable of recognizing xenogeneic porcine aortic endothelial cells in a calcium-dependent manner. Transcriptome-wide analysis to identify the differentially expressed transcripts in these cells revealed that the most likely candidate of the Galα1,3-Gal-independent recognition moiety is the tetraspanin CD82. Abs to CD82 inhibited the calcium response and the subsequent activation invoked by xenogeneic encounter. Our data identify CD82 on innate immune cells as a major "xenogenicity sensor" and open new avenues of intervention to making xenotransplantation a clinical reality.

Global assessment of GU-rich regulatory content and function in the human transcriptome.

Unlike AU-rich elements (AREs) that are largely present in the 3'UTRs of many unstable mammalian mRNAs, the function and abundance of GU-rich elements (GREs) are poorly understood. We performed a genome-wide analysis and found that at least 5% of human genes contain GREs in their 3'UTRs with functional over-representation in genes involved in transcription, nucleic acid metabolism, developmental processes, and neurogenesis. GREs have similar sequence clustering patterns with AREs such as overlapping GUUUG pentamers and enrichment in 3'UTRs. Functional analysis using T-cell mRNA expression microarray data confirms correlation with mRNA destabilization. Reporter assays show that compared to AREs the ability of GREs to destabilize mRNA is modest and does not increase with the increasing number of overlapping pentamers. Naturally occurring GREs within U-rich contexts were more potent in destabilizing GFP reporter mRNAs than synthetic GREs with perfectly overlapping pentamers. Overall, we find that GREs bear a resemblance to AREs in sequence patterns but they regulate a different repertoire of genes and have different dynamics of mRNA decay. A dedicated resource on all GRE-containing genes of the human, mouse and rat genomes can be found at brp.kfshrc.edu.sa/GredOrg.

AU-rich elements regulate Drosophila gene expression.

In mammals, AU-rich elements (AREs) are critical regulators of mRNA turnover. They recruit ARE-binding proteins that inhibit or stimulate rapid mRNA degradation in response to stress or developmental cues. Using a bioinformatics approach, we have identified AREs in Drosophila melanogaster 3' untranslated regions and validated their cross-species conservation in distant Drosophila genomes. We have generated a Drosophila ARE database (D-ARED) and established that about 16% of D. melanogaster genes contain the mammalian ARE signature, an AUUUA pentamer in an A/U-rich context. Using candidate ARE genes, we show that Drosophila AREs stimulate reporter mRNA decay in cultured cells and in the physiological context of the immune response in D. melanogaster. In addition, we found that the conserved ARE-binding protein Tis11 regulates temporal gene expression through ARE-mediated decay (AMD) in D. melanogaster. Our work reveals that AREs are conserved and functional cis regulators of mRNA decay in Drosophila and highlights this organism as a novel model system to unravel in vivo the contribution of AMD to various processes.

ARED Organism: expansion of ARED reveals AU-rich element cluster variations between human and mouse.

ARED Organism represents the expansion of the adenylate uridylate (AU)-rich element (ARE)-containing human mRNA database into the transcriptomes of mouse and rat. As a result, we performed quantitative assessment of ARE conservation in human, mouse and rat transcripts. We found that a significant proportion ( approximately 25%) of human genes differ in their ARE patterns from mouse and rat transcripts. ARED-Integrated, another updated and expanded version of ARED, is a compilation of ARED versions 1.0 to 3.0 and updated version 4.0 that is devoted to human mRNAs. Thus, ARED-Integrated and ARED-Organism databases, both publicly available at http://brp.kfshrc.edu.sa/ARED, offer scientists a comprehensive view of AREs in the human transcriptome and the ability to study the comparative genomics of AREs in model organisms. This ultimately will help in inferring the biological consequences of ARE variation in these key animal models as opposed to humans, particularly, in relationships to the role of RNA stability in disease.

Integrated analysis of experimental data sets reveals many novel promoters in 1% of the human genome.

The regulation of transcriptional initiation in the human genome is a critical component of global gene regulation, but a complete catalog of human promoters currently does not exist. In order to identify regulatory regions, we developed four computational methods to integrate 129 sets of ENCODE-wide chromatin immunoprecipitation data. They collectively predicted 1393 regions. Roughly 47% of the regions were unique to one method, as each method makes different assumptions about the data. Overall, predicted regions tend to localize to highly conserved, DNase I hypersensitive, and actively transcribed regions in the genome. Interestingly, a significant portion of the regions overlaps with annotated 3'-UTRs, suggesting that some of them might regulate anti-sense transcription. The majority of the predicted regions are >2 kb away from the 5'-ends of previously annotated human cDNAs and hence are novel. These novel regions may regulate unannotated transcripts or may represent new alternative transcription start sites of known genes. We tested 163 such regions for promoter activity in four cell lines using transient transfection assays, and 25% of them showed transcriptional activity above background in at least one cell line. We also performed 5'-RACE experiments on 62 novel regions, and 76% of the regions were associated with the 5'-ends of at least two RACE products. Our results suggest that there are at least 35% more functional promoters in the human genome than currently annotated.

Analysis of overrepresented motifs in human core promoters reveals dual regulatory roles of YY1.

A set of 723 high-quality human core promoter sequences were compiled and analyzed for overrepresented motifs. Beside the two well-characterized core promoter motifs (TATA and Inr), several known motifs (YY1, Sp1, NRF-1, NRF-2, CAAT, and CREB) and one potentially new motif (motif8) were found. Interestingly, YY1 and motif8 mostly reside immediately downstream from the TSS. In particular, the YY1 motif occurs primarily in genes with 5'-UTRs shorter than 40 base pairs (bp) and its locations coincide with the translation start site. We verified that the YY1 motif is bound by YY1 in vitro. We then performed detailed analysis on YY1 chromatin immunoprecipitation data with a whole-genome human promoter microarray (ChIP-chip) and revealed that the thus identified promoters in HeLa cells were highly enriched with the YY1 motif. Moreover, the motif overlapped with the translation start sites on the plus strand of a group of genes, many with short 5'-UTRs, and with the transcription start sites on the minus strand of another distinct group of genes; together, the two groups of genes accounted for the majority of the YY1-bound promoters in the ChIP-chip data. Furthermore, the first group of genes was highly enriched in the functional categories of ribosomal proteins and nuclear-encoded mitochondria proteins. We suggest that the YY1 motif plays a dual role in both transcription and translation initiation of these genes. We also discuss the evolutionary advantages of housing a transcriptional element inside the transcript in terms of the migration of these genes in the human genome.

PromoSer: improvements to the algorithm, visualization and accessibility.

PromoSer is a web service that provides an easy and efficient approach to the batch retrieval of a large number of proximal promoters. Since its introduction last year, it has undergone continued development and expansion. At the core, there have been improvements in the filtering of the raw mRNA/EST sequences upon which all predictions are built, improvements in the alignments clustering and transcription start site prediction algorithms, and improvements in the backing database for increased performance. At the user interface level, improvements include enhanced functionality and user options, better integration with other resources on the web and a new visualization tool. PromoSer now also supports queries using a SOAP-based interface and XML-based responses. The service is publicly available at http://biowulf.bu.edu/zlab/PromoSer.

SeqVISTA: a new module of integrated computational tools for studying transcriptional regulation.

Transcriptional regulation is one of the most basic regulatory mechanisms in the cell. The accumulation of multiple metazoan genome sequences and the advent of high-throughput experimental techniques have motivated the development of a large number of bioinformatics methods for the detection of regulatory motifs. The regulatory process is extremely complex and individual computational algorithms typically have very limited success in genome-scale studies. Here, we argue the importance of integrating multiple computational algorithms and present an infrastructure that integrates eight web services covering key areas of transcriptional regulation. We have adopted the client-side integration technology and built a consistent input and output environment with a versatile visualization tool named SeqVISTA. The infrastructure will allow for easy integration of gene regulation analysis software that is scattered over the Internet. It will also enable bench biologists to perform an arsenal of analysis using cutting-edge methods in a familiar environment and bioinformatics researchers to focus on developing new algorithms without the need to invest substantial effort on complex pre- or post-processors. SeqVISTA is freely available to academic users and can be launched online at http://zlab.bu.edu/SeqVISTA/web.jnlp, provided that Java Web Start has been installed. In addition, a stand-alone version of the program can be downloaded and run locally. It can be obtained at http://zlab.bu.edu/SeqVISTA.

Site2genome: locating short DNA sequences in whole genomes.

SUMMARY: Many biological papers describe short, functional DNA sites without specifying their exact positions in the genome. We have developed a Web server that automates the tedious task of locating such sites in eukaryotic genomes, thus giving access to the context of rich annotations that are increasingly available for genome sequences. AVAILABILITY: http://zlab.bu.edu/site2genome/

PromoSer: A large-scale mammalian promoter and transcription start site identification service.

Proximal promoters have a major impact on transcriptional regulation. Studies of the sequence-based nature of this regulation usually require collection of proximal promoter sequences for large sets of co-regulated genes. We report a newly implemented web service that facilitates extraction of user specified regions around the transcription start site of all annotated human, mouse or rat genes. The transcription start sites have been identified computationally by considering alignments of a large number of partial and full-length mRNA sequences to genomic DNA, with provision for alternative promoters. The service is publicly available at http://biowulf.bu.edu/zlab/PromoSer/.

Learning curves for Gaussian process regression: approximations and bounds.

We consider the problem of calculating learning curves (i.e., average generalization performance) of gaussian processes used for regression. On the basis of a simple expression for the generalization error, in terms of the eigenvalue decomposition of the covariance function, we derive a number of approximation schemes. We identify where these become exact and compare with existing bounds on learning curves; the new approximations, which can be used for any input space dimension, generally get substantially closer to the truth. We also study possible improvements to our approximations. Finally, we use a simple exactly solvable learning scenario to show that there are limits of principle on the quality of approximations and bounds expressible solely in terms of the eigenvalue spectrum of the covariance function.